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1. Introduction 

New approximate methods of nuclear structure theory are usually examined by applying 
them to simple exactly soluble models in order to gain some insights on a range of their 
validity. One of the widely used models is the two level schematic shell model which 
possess the 577(2) symmetry and is often called SU(2) or Lipkin - Meshkov - Glick 
(LMG) model [1]. Numerous applications of the LMG- model can be found in [2]. 

During the last years the model has been used many times to justify approximate 
methods of the many-body theory at finite temperature [3-9]. These methods are espe- 
cially interesting in view of current intensive studies of hot nuclear systems. Previous 
works on the LMG- model at finite temperature [3,4] have focused on boson expansion 
methods and symmetry breaking in hot nuclei. The mixed state representation has been 
formulated and then applied to the LMG- model in refs.[5-8]. In particular, the ther- 
mal Hartree - Fock approximation as well as the thermal random phase approximation 
(TRPA) were studied within the approach [8]. The thermal Hartree - Fock approximation 
and the static path approximation were analyzed within the model in ref. [9]. 

A new approximate method has been recently proposed [10] to describe collective 
excitations in hot finite Fermi systems. This method, the so-called renormalized TRPA 
(TRRPA), is an extension of the renormalized RPA of Ken-ji Hara [11] and D. Rowe 
[12] to finite temperatures. Within TRRPA vibrational excitations are supposed to be 
harmonic like in TRPA but a temperature - dependent ground state is treated in a more 
consistent manner. Namely, a finite number of thermal quasiparticles are presented in 
this ground state. 

In the present paper, we investigate the accuracy of the thermal renormalized random 
phase approximation by comparing it with the exact calculations for the grand canonical 
ensemble for the LMG- model. Moreover, a comparison with the thermal mean field 
approximation (TMFA) and TRPA is also done. 

2. The LMG- model and the grand canonical ensem- 
ble calculations 

The following version of the LMG- model is used: N fermions are distributed over two 
levels, each having a degeneracy f2. The distance between the levels is e, the coupling 
constant V does not depend on any quantum number. At T = and V = the lower 
level is full, the upper - empty, i.e. N = £1. The model Hamiltonian has the form 




(1) 



where the operators of quasispin J and its projections J + , J_, J z are defined as follows: 

J 2 = I (j + J+ + J-J-) + Jt , 
1 n n 



p=l p=l 

Here a J and a ip are particle creation and annihilation operators on the lower {i = 1) or 
the upper (i = 2) level. 

The operators J± and J z form SU(2) algebra, and the quasispin operator commutes 
with H . So the Hamiltonian matrix breaks up into submatrices Oj of dimension 2 J + 1. 
The Hamiltonian can be diagonalized in each of these subspaces independently. The 
corresponding eigenvalues are denoted by E(, E^-, ...-E/j+i- 

To calculate the grand canonical partition function, one needs the eigenvalues E'l and 
the degeneracies of irreducible quasispin representations Oj for different particle numbers 
from the range < iV < 2Q. The latter have been determined in ref. [7], and here we 
use this result. The whole number of the ensemble states, i.e., the whole number of the 
eigenstates of the LMG- systems formed by two Q- degenerated levels with a number of 
particles varying from 1 to 2Q is equal to 2 2Q . Any state of the ensemble can be written 
in the following form: 

\giPl,92P2, -gnPn) = < ipi O+ P2 , •••O+pJO) , a 9p |0) = , 

The indices gi,pi,n have the following meanings: 

gi E {1,2}, Pi e % e {l,...n}, n e {l,...2fi}, 

i.e., g marks the lower and the upper levels, p - sublevels, i is an index of a particle and n 
is the particle number in the particular LMG- system from the grand canonical ensemble. 
If n — 0, \g\Pi, <?2P2, ■■■9nPn) = |0). A particular distribution of the given number of 
particles over two degenerate levels can be characterized by numbers v\ and v^. v\ is a 
number of sublevels which are occupied by particles for both the lower and upper levels; 
v<i is a number of sublevels which are occupied for neither the lower nor the upper level. 
The quasispin J of the state is determined by the distribution of the rest of particles over 
It sublevels where 2r = Q — v\ — 1*2- The number 2(r + v\) is equal to the number of 
particles. We denote by r piiP2) ... P2T+ ^ i the subspace of states with v\ occupied and v 2 empty 
sublevels. Its dimension is 2 2t . There exist f2!/(2r)!z/i!z/2! distinct subspaces r piiP2) ...p 2T+i/i 
for fixed r and v\. Each of them may be decomposed into irreducible subspaces with 
fixed quasispin values T (appearing once), T _i (appearing g\ times), O r -2 (appearing 
g\ times), ... , Q T -k (appearing g\ times),... , 6 T _[ T ] (appearing g^ times). Here 

(2r)! (2r)! 
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k\{2r-k)\ (k - l)\(2r - k + iy: 



and [t] = r, if r is integer, [r] = r — 1/2 if r is half-integer. 
Then, the exact grand partition function is 



rvif2 V / 1 ^ fc rra 



^- fc -2(r + ^)A 
T 



The expressions for average energy, quasispin z-projection and the total fermion number 
are 
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3. Thermo field dynamics: basic elements 

To be more understandable while describing approximate methods, we briefly recapitulate 
the formalism of thermo field dynamics (TFD) (see, [3, 13-15]). 

The extension of quantum field theory to finite temperature requires the field degrees 
of freedom to be doubled. In TFD, this doubling is achieved by introducing an additional 
tilde space. A tilde conjugate operator A is assigned to an operator A (acting in ordinary 
space) through the tilde conjugation rules 

(AB) = AB ; (aA^+bB) = a* A + b*B , 

where A and B represent ordinary operators and a and b are c-numbers. The asterisk 
denotes the complex conjugation. The tilde operation commutes with hermitian conju- 
gation and any tilde and non-tilde operators are assumed to commute or anticommute 
with each other. A double application of tilde operation changes a sign of a fermionic 
operator and saves it for a bosonic one. The whole Hilbert space of a heated system is a 
direct product of ordinary and tilde spaces. A formal quantity playing a central role in 
the present discussion is the so-called thermal Hamiltonian: 

h = h-h 

The operator 7i serves to translate temperature dependent wave functions along the time 
axis. It means that an "excitation spectrum" of a hot system (or, in other words, a set 
of energies corresponding to the thermal equilibrium states) should be obtained by the 
diagonalization of 7i. 



The temperature-dependent vacuum |\I/ (T)) is the eigenvector of TC with eigenvalue 



H\*o(T)) = 
If one determines the thermal vacuum state as 

\MT)) = , -- = — = E«p(-^)1") ® l»> 

where |n) and |n) are eigenvalues, eigenvectors and their tilde counterparts of the 
Hamiltonian respectively, the expectation value (\E'o(^)|0| , I / o(^)) will exactly corre- 
spond to the grand canonical ensemble average <C O ^> of a given observable O. 

In practice, it is impossible to find the exact thermal vacuum for a full Hamiltonian 
of a many-body system. In setting up approximate schemes, the usual starting point is 
the thermal mean-field approximation. In this case, the thermal vacuum \^ (T)) is an 
eigenvector of the uncorrelated thermal Hamiltonian 

Hmf\0(T)) = (H MF -H MF )\0(T))=Y^ei(afa i -ata i )\0(T)) = 0. (2) 

i 

The solutions of eq. (2) define the vacuum |0(T)) for so-called thermal quasiparticles (3, (3: 

fa = x&i + Didf 

A|0(T)) = A|0(T)) = 0, 

where the coefficients Xj, denote the thermal Fermi occupation probabilities of the states 
a i"|0) (|0) is a vacuum for a,) 



x% = v 1 - ft , Vi = v L 
1 



1 + exp(^/T) 

Sometimes the {x, y} transformation is called the thermal Bogoliubov transformation. It 
is a unitary transformation and thus conserves the commutation relations. 



4. Approximate methods 

Now we apply the TFD formalism to the LMG- model and derive the corresponding 
equations of TRRPA. A more general formulation of the thermal renormalized random 
phase approximation can be found in refs. [10,16,17]. Moreover, within the Hartree - Fock 
method, depending on the value of the coupling constant V two different phases of the 



LMG- system exist - a normal phase and a deformed one. The present consideration is 
restricted to a normal phase. So we do not take into account the mean field rearrangement 
which occurs if the value of the effective coupling constant x = V(N — becomes more 
than unity. 

The model thermal Hamiltonian 7i = H — H, where H has the form (1), has to be 
written in terms of the thermal quasiparticle operators. The first item in (1) conserves 
the diagonal form. The interaction operator becomes more complicated. For further 
studies we need only that part of 7i which consists of the terms with even numbers of 
both creation and annihilation thermal quasiparticle operators. Namely, 

v{h - / 2 ) 



H' = e(B-B)- 



(A +2 + A 2 ) - (i+ 2 + A 2 ) 



(3) 
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where 

B = \ E K&p - PtAp) > A+ = E W P , ^ + = E W 
p=i p=i p=i 

The following exact commutation rules are valid for the thermal biquasiparticle oper- 
ators A, A + , A and A + : 

a, a + ] = n - f: M P hv - E tiv^v , [a, a + ] = n - jr p+p lp - M P K ■ (4) 



P =i 



P =i 



P =i 



P =i 



All other commutators between the operators A, A + , A and A + vanish. 

By the use of the Wick theorem one can approximate [10,16] the r.h.s. of (4) by 
c-numbers neglecting an influence of the pairs of normal ordered operators : f3 + f3 : and 
: (3 + (3 :. Namely, 



A, A* 



A,A + =N(l- Pl -p 2 ) = N(l-2p) . 



(5) 



Here pi are the numbers of thermal quasiparticles in the temperature - dependent ground 
state \^>o(T)) that will be defined later on. That is 

P* = jf(MT)\N? \MT)) = ^(MT)\N?\MT)) 

where TVf is the operator of the number of thermal quasiparticles iVf = X)p=i Pipfrp ■ 

The thermal Hamiltonian (3) can be diagonalized in the space of two one-phonon 
states constructed as bilinear forms of the thermal biquasiparticle operators: 



g+|M/ o (T)) = (^A + -0 1 A)|M/ o (T)) 

Qt\Vo(T)) = U 2 A + - faA) |^o(T)) . 



(6) 



Now we define the wave function of the temperature - dependent ground state \^q(T)) as 
the thermal phonon vacuum, i.e. Qi^\^q{T)) = 0. 



The states (6) have to be orthonormal. Thus, taking account of eq. (5) the following 
constraints on the amplitudes ip and <fi are derived 



- 2 = [N(l - 2P)]- 1 , 1 = 1,2. 

The system of equations for fy, (pi and the phonon frequencies iO{ is easily obtained by the 
equation of motion method. It appears that only a positive value of lo\ and a negative value 
of u 2 is allowed under a requirement that the wave functions Qi\^/ (T)) and Q 2 \^o(T)) 
are vectors of the Hilbert space. The eigenvalue - eigenvector problem has the following 
solution: 

Ul =u> = V / ^-V 2 (/2-/i) 2 (l-2p) 2 (iV-l) 2 , 

, 2 e + 0J , 2 e — UJ 



2Nu{l - 2p) ' ri 2Nu(l - 2p) ' 

u 2 = -w, i>l = i>l, 4>l = 4>\- 

One more equation has to be added to the above system - the equation for p. To 
evaluate this equation we need an expression for the thermal phonon vacuum state. The 
latter can be derived from the thermal quasiparticle vacuum state |0(T)) by a unitary 
transformation 

|*„(T)> = VRe s \0(T)) , S = 2(1 ^ 2p) ^ {A + A + + A+A+) . 
By the use of standard techniques of the operator calculus [11] we get 

It is interesting to note that in the thermodynamic limit, i.e. at N — > oo, p vanishes and 
the TRRPA equations are reduced to the TRPA ones. 

Let us display the expressions for the average energy, the average quasispin z-projection 
and the variance of the particle number 

\ti )trrpa — ~ (1 - 2p) H x 



2 2cu 2(/ 2 -/0 



\Jz)trrpa — [i- - 2p) 



AN TRRPA = sjN{l-2p) [/i(l-/i) + / 2 (l-/ 2 )] 



The above expectation values were taken over the TRRPA vacuum state. 

It seems appropriate to give expressions for the same quantities within other approxi- 
mations - TRPA and TMFA. The TRPA expressions are obtained from the TRRPA ones 



by putting p = 0. In this case, the commutator [A, A + ] is equal to N and the expectation 
values are taken over the TRPA vacuum 

tfrn Ne{f 2 -h) e^-uP (h-fi? + l 

\-n /TRPA = ~ 1 ~ x 



(Jz) 



TRPA 



2 2cu 2(/ 2 - h) 

N(f 2 - h) 



AN TRPA = y/N [/i(l-/i) + / 2 (l-/ 2 )] • 



Within TMFA the interaction between particles is omitted and one has to evaluate the 
expectation values over the thermal quasiparticle vacuum |0(T)). The TMFA expressions 
for (J z )tmfa and ANtmfa are the same as in TRPA. For (H)tmfa one gets 

( M Ne(f 2 - h) 

\rt)TMFA — ~ • 



5. Results and discussion 

The numerical calculations are done for the LMG- system with N = £1 =10 particles and 
e—2. The results are displayed in Figs. 1-6. 

Firstly, we discuss a dependence of some characteristics of the system on the effective 
coupling constant x- The energy of the collective state u as a function of x at T = and 
0.25e is displayed in Fig. I. Besides the results of the TRPA and TRRPA calculations 
the exact solution at T =0 is also shown. As it should be, with increasing x the energy uo 
goes down. Within TRPA uo vanishes at x — 1- This collapse does not take place for the 
exact solution as well as for the TRRPA result. This last feature of the RRPA solution 
is well known in the case of a cold nucleus and is actively used in some recent nuclear 
structure calculations [18]. As it has been demonstrated for the first time in ref.[10], the 
same is valid at T ^ 0. In the present version of the LMG- model heating effectively 
weakens the interaction of particles (the effective coupling constant x is multiplied by a 
factor of fi — f 2 < 1) and the TRPA collapse occurs at larger x~ values. In TRPA when 
X Xcoiiapse (H)trpa ~ —uj^ 1 — > — oo. It is not the case for TRRPA (see Fig. 2). The 
value (H)trrpa goes down much slower and remains even greater than the exact value 
(H)gce- At large values of x the strong difference between {H)trrpa and {H)gce is due 
to neglecting the mean field rearranging. 

In Figs. 3-5, the average energy of the system, the average J z value and the variance of 
the particle number as functions of T are displayed (x =0,95). The noticeable difference 
between the exact and the approximate values is only at moderate T < 0, 3e. Here 
TRRPA works evidently better than TRPA and TMFA. The absolute values {H) TRP a 



and (J z )trpa are greater than (H) GC e and (J z )gce, respectively. At the same time, 
\(H) T rrpa\ < \{H)gce\- The relation | ( J z ) T rrpa\ < \{J z )gce\ is valid only at T < 0, 8e 
but then \(J z )trrpa\ appears to be slightly greater than \(J z )gce\- At T > 0,5^ the 
differences between the exact and the approximate results is negligible. The difference 
between exact and approximate values of the particle number variance is only 2-3%, i.e. 
even less than for other variables. Decreasing in the difference with raising up T is a 
result of effective weakening of the interaction. 

The value AN/N as a function of N is shown in Fig. 6. It decreases slowly when T 
increases, and its typical value at N = 10-30 is around 10%. The approximate methods 
disturb AN only slightly. The difference between different approximations is of minor 
importance although formally TRRPA seems to be better. 

6. Summary 

Taking the Lipkin - Meshkov - Glick model as an example we have studied a validity 
of some approximate methods of many-body theory at finite temperature. The average 
energy, the average quasispin z-projection and the particle number variance as functions 
of temperature and particle number have been calculated in different approximations 
as well as exactly with the grand canonical partition function. On the whole, TRRPA 
gives better results than other approximations. Its advantages are especially evident at 
moderate temperatures T < 0, 5e. With increasing T and N, results of approximate 
methods improve rapidly and at T e the difference between exact and approximate 
results is invisible. 

In the present paper, we have studied only the case with not too strong particle 
interaction (x < 1). Investigations of the deformed phase of the LMG- model are in 
progress. 
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Figure captions 



Fig. 1 The energy of the lowest excited state in the LMG- model as a function of the 
effective coupling constant x & t T = and 0,25e. The exact results - open circles; 
the TRPA (RPA) results - dashed lines; the TRRPA (RRPA) results - solid lines. 

Fig. 2 The average energy of the LMG- system (H) as a function of the effective coupling 
constant X- The exact results for the grand canonical ensemble - open circles; the 
TRPA results - dashed line; the TRRPA results - solid line. 

Fig. 3 The average energy (if) as a function of temperature T. The exact results for 
the grand canonical ensemble - open circles; the TRPA results - dashed line; the 
TRRPA results - solid line. 

Fig. 4 The average value of the quasispin projection (J z ) as a function of temperature T. 
For notation see Fig. 3. 

Fig. 5 The particle number variance AN as a function of temperature T. For notation see 
Fig. 3. 

Fig. 6 The dependence of AN/N on a particle number N. For notation see Fig. 3. 
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